The large degeneracy of excited hadrons and quark models 
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The pattern of a large approximate degeneracy of the excited hadron spectra (larger than the 
chiral restoration degeneracy) is present in the recent experimental report of Bugg. Here we try to 
model this degeneracy with state of the art quark models. We review how the Coulomb Gauge chiral 
invariant and confining Bethe-Salpeter equation simplifies in the case of very excited quark-antiquark 
mesons, including angular or radial excitations, to a Salpeter equation with an ultrarelativistic 
kinetic energy with the spin-independent part of the potential. The resulting meson spectrum 
is solved, and the excited chiral restoration is recovered, for all mesons with J > 0. Applying the 
ultrarelativistic simplification to a linear equal-time potential, linear Regge trajectories are obtained, 
for both angular and radial excitations. The spectrum is also compared with the semi-classical Bohr- 
Sommerfeld quantization relation. However the excited angular and radial spectra do not coincide 
exactly. We then search, with the classical Bertrand theorem, for central potentials producing always 
classical closed orbits with the ultrarelativistic kinetic energy. We find that no such potential exists, 
and this implies that no exact larger degeneracy can be obtained in our equal-time framework, with a 
single principal quantum number comparable to the non-relativistic Coulomb or harmonic oscillator 
potentials. Nevertheless we find plausible that the large experimental approximate degeneracy will 
be modelled in the future by quark models beyond the present state of the art. 



I. INTRODUCTION 

In the recent report of Bugg 0, H, H, [3 a large de- 
generacy emerges from the spectra of the angularly and 
radially excited resonances produced in pp annihilation 
by the Crystal Ball collaboration at LEAR in CERN 
This degeneracy may be the third remarkable pattern of 
the excited spectra of hadronic resonances. 

A long time ago, Chew an Fautschi remarked the exis- 
tence of linear Regge trajectories [6[ for angularly excited 
mesons, and baryons 0, A similar linear aligning of 
excited resonances was also reported for radial excita- 
tions 0. 

Recently, Glozman et al. [!, [ll[ have been system- 
atically researching the degeneracy of chiral partners in 
excited resonances, both in models and in lattice QCD. 
Le Yaouanc et al. [12] and PB et al. (15 , Il4j developed a 
model of spontaneous symmetry breaking. Viewed retro- 
spectively, this already included the degeneracy of chiral 
partners in the limit of high radial or angular excitations 
[15| , both in light-light and in heavy-light hadrons. This 
earlier work was based on the Bogoliubov transforma- 
tion and used the Coulomb gauge or the local coordinate 
gauge QCD truncations. Here the formalism will be ex- 
tended to include spin-dependent interactions explicitly. 

The recent data on highly excited mesons, observed by 
LEAR at CERN [H,0,[5[ and of excited baryons observed 
by the Crystal Barrel collaboration at ELSA stress 
the interest of possible patterns in the excited hadronic 
resonance spectra. While we still have to wait for new ex- 
periments focused on excited mesonic resonances to con- 
firm the report of Bugg, say in PANDA, GLUEX, or 
BESIII, it is important to research theoretical models of 
the excited hadrons. 

Here we address the question, is it possible to build 
an equal-time quark model, with linear trajectories, with 
excited chiral symmetry, and, also, with a principal quan- 



tum number? We adopt the framework of the Coulomb 
gauge confinement, of the mass gap equation, and of 
the equal-time Bethe-Salpeter equation. In Section II 
we review and expand earlier work on chiral symmetry 
breaking and mesonic boundstates, and show in detail 
the equations for the simplest potential. We also show 
how for excited states (and for J > 0) the Bethe-Salpeter 
equation simplifies to a Schrodinger-like Salpeter equa- 
tion, with ultrarelativistic (massless) kinetic energies and 
a chiral symmetric equal-time potential. In Section III 
we solve the equation with the method of the double 
diagonalization of the equal time hamiltonian and show 
how linear equal-time potentials and massless quarks pro- 
duce linear Regge trajectories, both for angular and ra- 
dial excitations. We also compare them with the Bohr- 
Sommerfeld semi-classical quantization. We then address 
in Section IV the large degeneracy, where both radial 
and angular excitations are degenerate. Extending to 
ultrarelativistic particles the techniques of the classical 
Bertrand theorem on closed orbits, we verify that no in- 
stantaneous 2-body potential may exactly produce the 
desired large degeneracy. Nevertleless we present plausi- 
ble solutions to this modelling problem in the conclusion, 
Section V. 



II. QUARK MASS GAP AND BOUNDSTATES 
IN EQUAL TIME 

We first review earlier work on chiral symmetry break- 
ing with equal-time confining quark-quark potentials, 
and show the example of the simplest possible model of 
this class of potentials, which continues to be explored 
[T3|. Importantly, the hamiltonian of this model can be 
approximately derived from QCD, 
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FIG. 1: The constituent quark masses m c (k), solutions of the 
mass gap equation, for different current quark masses mo. 
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up to the first cumulant order, of two gluons 

|20|. In the modified coordinate gauge the cumulant is, 

g 2 (A a ^x)AUy)) * -^<W>3,o [^ 3 (x - y) 2 - U] (2) 

and this is a simple density-density harmonic effective 
confining interaction, ttiq is the current mass of the 
quark. The infrared constant U confines the quarks but 
the meson spectrum is completely insensitive to it. The 
important parameter is the potential strength Kq, the 
only physical scale in the interaction, and all results can 
be expressed in units of K$. A reasonable fit of the 
hadron spectra is achieved with Kq ~ 0.3 ± 0.05 GeV. 
The relativistic invariant Dirac-Feynman propagators 
can be decomposed in the quark and antiquark 
Bethe-Goldstone propagators [2lj . used in the formalism 
of non- relativistic quark models, 
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where S = sin(y) 



C = cos(<^>) 



and tp is a chiral angle. In the non condensed vac- 
uum, ip is equal to arctan ^ . In the physical vacuum, 
the constituent quark mass m c {k), or the chiral angle 
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is a variational function which is 



determined by the mass gap equation. We anticipate 
examples of solutions, for different light current quark 
masses too, depicted in Fig. Q] 

There are three equivalent methods to derive the mass 
gap equation for the true and stable vacuum, where con- 
stituent quarks acquire the constituent mass [l7j . One 
method consists in assuming a quark-antiquark 3 P con- 
densed vacuum, and in minimizing the vacuum energy 
density. A second method consists in rotating the quark 
and antiquark fields with a Bogoliubov-Valatin canonical 
transformation to diagonalize the terms in the hamilto- 
nian with two quark or antiquark second quantized fields. 
A third method consists in solving the Schwinger-Dyson 
equations for the propagators. Any of these methods lead 
to the same mass gap equation and quark dispersion re- 
lation. Here we replace the propagator of eq. (J3]) in the 
Schwinger-Dyson equation, 
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where, with the simple density-density harmonic interac- 
tion [12j, the integral of the potential is a laplacian and 
the mass gap equation and the quark energy are finally, 
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Numerically, this equation is a non-linear ordinary differ- 
ential equation. It can be solved with the Runge-Kutta 
and shooting method. Examples of solutions for the cur- 
rent quark mass m c (k) = ktantp, for different current 
quark masses mo, are depicted in Fig. [T] 

The Salpeter-RPA equations for a meson (a colour sin- 
glet quark-antiquark bound state) can be derived from 
the Lippman-Schwinger equations for a quark and an an- 
tiquark, or replacing the propagator of eq. in the 
Bethe-Salpeter equation. In either way, one gets [2l[ 
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TABLE I: The positive and negative energy spin- independent, 
spin-spin, spin-orbit and tensor potentials, computed exactly 
in the framework of the simple density-density harmonic 
model of eq. ©. <p'(k), C(k) and Q(k) = 1 - S(k) are all 
functions of the constituent quark(antiquark) mass. 
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where fci 



and P is the total 



ci = k + , k 2 = k - 
momentum of the meson. 

The Salpeter-RPA equations of PB et al. T3] and 
of Llanes-Estrada et al. [H| are obtained deriving the 
equation for the positive energy wavefunction <f) + and for 
the negative energy wavefunction <j>~ . Solving for one 
gets the Salpeter equations of Le Yaouanc et al. 
This results in four potentials V a @ respectively coupling 
v a = r<fi a to v" , in the boundstate Salpeter equation, 
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The relativistic equal time equations have the double of 
coupled equations than the Schrodinger equation. The 
negative energy component v~ is smaller than the pos- 
itive energy component by a factor of the order of 1 /M 
in units ok Kq = 1. Thus when M is large, and this 
is the case for most excited mesons, the negative energy 
components can be neglected and the Salpeter equation 
simplifies to a Schrodinger equation. 

Importantly, the potentials V ++ = V and V + ~ = 
V~ + include the usual spin-tensor potentials [23], pro- 
duced by the Pauli a matrices in the spinors of eq. ([3]). 
They are detailed explicitly in Table [JJ Because we are 
interested in highly excited states, where both (r) and 
(fc) are large, we consider the limit where ^ — ► 0. This 
implies that the potentials, used in Table [I] f'(k) — > 0, 
C(k) -> 1 and Q(k) = 1 - S(k) -> 1. Then using the text- 
book matrix elements of the spin-tensor potentials, the 
boundstate Salpeter equation decouples in two different 
equations depending only on J and not explicitly on L 
or S. Without the chiral degeneracy there would be four 
different resonances for each j, one with s = and j = I 
and three with s = 1 and j =1 — 1, I, l + l. With the chi- 
ral degeneracy we get only two different equations, one 
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TABLE II: Masses of the first angular and radial excitations 
of the different light-light tachyons and mesons in the chiral 
limit of m << k. Each column includes both positive and 
negative parity degenerate states, except for the pseudoscalar 
and scalar tachyonic states. The meson masses are separated 
in two different families with the same J because two different 
Salpeter equations ([5| and @ exist for each J. 
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and another for j ' > 1 with, 
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Thus states with different I and equal j, i.e. with different 
parity are degenerate, and chiral symmetry is restored. 
This chiral degeneracy applies to all angular momenta, 
except for j = 0. In Table [TT] we show the masses of 
the different light-light mesonic solutions of eqs. (jS} and 
([5]), including tachyons, corresponding to the limit where 
m c << k. 

The j = case is a subtle one, because neglecting the 
mass m c is equivalent to consider the chiral limit, and 
this is equivalent to change from the physical vacuum 
to the chiral invariant vacuum which is a false, unstable 
vacuum. A detailed inspection shows that the different 
potentials —j^z , 2fc , -p- are bound from below and pos- 
itive definite in the sense that all their eigenvalues are 
positive. However — p- is unbound from below. It turns 
out that for j = all the solutions of eq. ([8|), including 
all radial excitations, are tachyons (24|, relevant for the 
structure of the chiral invariant false vacuum of QCD 
[25[, corresponding to a different solution of the mass 
gap equation m c = 0. Even when a very small regular- 
izing quark mass m c is assumed, constant for simplicity, 
the tachyons persist. This is confirmed by the numerical 
solutions of the regularized Salpeter equation, shown in 
Table [TT] . Technically, in the j = case, it is necessary 
to rescale the momentum and mass, 
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where any finite solution M' in fact corresponds to a large 
mass M — M'/rri^, and where a wave- function with a 
finite k' corresponds to a wave-function with small mo- 
mentum k = k' m c . A long time ago Le Yaouanc et 
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al. [T3| showed that in the chiral limit the pseudoscalar 
and the scalar possess tachyonic solutions. Very rece ntly 
PB showed that this number of tachyons is infinite [241 ]. 
Only with a finite m c quark mass, do the scalar and pseu- 
doscalar mesons have positive masses. But then the ex- 
cited scalar and pseudoscalar states are not degenerate, 
in contradistinction with the chiral degeneracy of the ex- 
cited mesons with j > 0. 

For excited mesons with j > 0, the spectrum in table 
HT1 is very well aproximated by the solutions of the pair 
of Schrodinger-like equations, 
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obtained when the negative energy component v 
neglected in eqs flSJ) and ©. 
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TABLE III: Masses of the light-light mesons, in dimension- 
less units of a = 1, computed with the ultrarelativistic equal 
time chiral degenerate Schrodinger equation (|13[) . The j — 
mesons are distant from the experimental spectrum, but chi- 
ral degeneracy is theoretically plausible for the very excited 



III. LINEAR REGGE TRAJECTORIES AND 
SEMI-CLASSICAL QUANTIZATION 

The simple quadratic model is now extended to a linear 
potential [22|, [26[, to get the linear Regge trajectories. 
It is well known, both from the quark modelling of the 
hadronic spectra, and from Lattice QCD, that the long 
range confining quark potential is linear. Notice that 
Szczepaniak et. al. 27], in the coulomb gauge, were able 
to derive from QCD the linear potential. Continuing with 
the limit where — > 0, we assume the radial equation, 
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dr 2 
Ev{r) , 



i(j + i) 
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where the negative energy components (vanishing with 
jj) were neglected, and where the terms similar to the 
terms -p- or ^ of eq. (fTTj) or eq. (fT2")) were also neglected. 
We only expect these approximations to be reasonable for 
highly excited states, and for j > because the number 
of excited false vacua remains infinite [28[ . 

Then we have to solve a Salpeter equation (or 
Schrodinger equation with ultrarelativistic kinetic en- 
ergy) except that the spherical angular momentum I is 
now replaced by the total angular momentum j, in the 
centrifugal barier. Eq. (fT3]) is solved with the method of 
the double diagonalization of the equal time hamiltonian. 
Using finite differences, first we diagonalize the bounded 
from below operator + ^{.t^ ■ Then we apply the 
square root. After returning to the original position space 
basis, we diagonalize the full hamiltonian. This provides 
us automatically with the full spectrum including the an- 
gularly and radially excited states. 

The results are shown in Table IIIII In Fig. [2] and in 
Fig. [3]we graphically demonstrate that the angular exci- 
tations and the radial excitations of this simple spectrum 



are disposed in linear Regge trajectories, 
j ~ a + aM 2 , 
j — Po + PM 2 ■ 



(14) 



This agrees qualitatively with the experimental spectrum 
[l|, where the linear Regge trajectories are also present. 

Interestingly, with the semi-classical Bohr-Sommerfeld 
quantization relation 



p dq ~ nh 



(15) 



and with the energy 2p c + a r of eq. (fT3"|) the linear 
trajectories can be derived. The linear trajectories for 
the angular excitations can be derived from the circular 
classical orbits, 

f 2 (2tt § p) ~ I h 
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where we also used the centripetal acceleration. The lin- 
ear trajectories corresponding to radial excitations can be 
derived from the linear classical orbits with L = 0, 
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Thus, in our units ofh = c = a = l we get for the Regge 
slopes respectively, 
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FIG. 2: We show the quasi-linear Regge trajectories, of j 
as a function of M 2 . Each line corresponds to a fixed radial 
n, increasing from left to right. The M are the masses of 
the light-light mesons, in dimensionless units of a — 1, com- 
puted with the ultrarelativistic equal time chiral degenerate 
Schrodinger equation (fl3j). 



FIG. 3: We show we show the quasi-linear Regge trajectories, 
of n as a function of M 2 . Each line corresponds to a fixed an- 
gular j, increasing from left to right. The M are the masses of 
the light-light mesons, in dimensionless units of a = 1, com- 
puted with the ultrarelativistic equal time chiral degenerate 
Schrodinger equation (fl3j). 



in excelent agreement with the slopes of Figs. [2] and [3J 
and similar to a recent Bethe-Salpeter calculation [lfj . 

Although the simplified model of eq. (fl3|) agrees qual- 
itatively with the experiment, both in the chiral degener- 
acy and in the linear Regge trajectories, there are quan- 
titative discrepancies. In particular, experimentally, the 
radial and agular slopes defined in cq. (fT4j) should be 
almost identical [l|, 

a ex p = 0.877 GeV~ 2 , 

/3exp = 0.855 GeV" 2 , (19) 

while our theoretical slopes, depicted in Figs. [2] and [3J 
and semi-classically computed in eq. (H"8"l) . differ by 7r/2. 
Moreover the meson masses in Table IIIII do not exactly 
comply with the large degeneracy emerging in the obser- 
vations of Bugg. 

IV. SEARCHING FOR A PRINCIPAL 
QUANTUM NUMBER WITH CLASSICAL 
CLOSED ORBITS 

To better model the large degeneracy, we enlarge our 
class of time- like central potentials, searching for the best 
possible potential of this class. Actually the model of 
Section III, with a linear potential only, is oversimplified. 
To fit the correct positive intercepts ao and @q of eq. 
(|14|) . it is standard in quark models to include in the 
potential a constant negative energy shift and a negative 
short range Coulomb potential. 

To exactly reproduce the large degeneracy, we ask for 
a spectrum with a principal quantum number, similar to 
the non-relativistic spectra of the Coulomb potential or 
of the harmonic oscillator potential, where the principal 
quantum numbers are respectively n + I + 1 and In + 



I + |. The difference here is that our kinetic energy is 
the ultrarelativistic one T = 2pc whereas in the non- 

2 

relativistic case T = 

Notice that in the non-relativistic case, classical closed 
orbits coincide with a quantum principal number. We 
can also address this problem searching for classi- 
cal closed orbits with an ultrarelativistic kinetic en- 
ergy. When all the classical orbits are closed then 
the Hermann- Jacobi-Laplace-Runge-Lenz vector is con- 
served, because this vector does not precess. In the 
Hamilton- Jacobi formalism, this vector commutes with 
the hamiltonian. The same commutation then also oc- 
curs in the quantum Schrodinger formalism, the formal- 
ism we are using now. Then a larger symmetry group, in- 
cluding the angular momentum and the Hermann- Jacobi- 
Laplace-Runge-Lenz vector exists. Finally this implies 
that a principal quantum number exists. 

Thus, rather than solving the ultrarelativistic 
Schrodinger equation for all the infinite possible central 
potentials, we prefer to extend the Bertand theorem tech- 
niques to the search of classical closed orbits with the ul- 
trarelativistic kinetic energy. For simplicity we consider a 
kinetic energy T = pc and a general potential V(r), used 
for a single particle in a central potential, comparable to 
our two-body problem in the centre of mass frame. 

Let us consider classical planar trajectories for an ul- 
trarelativistic quark with the speed of light, and with 
momentum, 

P 

p = -v 

c 

= -re r + ?rfl? 9 , (20) 
c c 

subject to a central force, 
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where the notation is obvious. In the plane we have two 
constants, the angular momentum L and the energy E, 



L = -r 2 9e ± . 

c 

E = pc + V(r) 



(22) 



Then Newton's law produces the equation for the radius 
as a function of time, 



E —V 



d 

dr 



(23) 



To study the condition to get closed orbits it is convenient 
to relace the variable time t by the polar coordinate 9 and 
the function r by it's inverse u = 1/r. Then the equation 
simplifies to, 



d9 2 ' 



J(u) 



J(u) 
E 



V(l/u) d 



du 



c 2 E 2 



V(l/u) 



(24) 



To get the classical case one only has to replace in the 
right hand side of eq. ([24|) the factor E ~ V c { 1 ^ u " > — > m. 
Thus the ultrarelativistic case has two independent con- 
stants E and cL while the classical case has only one 
constant m/L. 

The theorem of Bertrand can be addressed considering 
a trajectory u close to the circular trajectory uq, 



u = u + rj 



(25) 



where we fix the angular momentum L and the trajectory 
is defined by the inverse radius u and derivative at 9 = 
and by the energy E. Defining, 



P 2 
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du 



J\u-- 
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we get, to leading order in r), 
d 2 r, 



d9 2 



+ P*T) = 



rj = hi cos (3(9 — Q ) 



(27) 



Thus for very small perturbations to the circular orbit, 
the condition for a closed orbit is that (3 is an integer 
number (if we want the orbit to close right after one turn) 
or rational (if we want it to close after a finite number 
of turns). Importantly this implies that (3 is constant, 
since it can't change continuously from one orbit with 
i*o to the next one. This restricts the class of possible 
potentials, 
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and the problem here is that, unlike in the non-relativistic 
limit, the equation still depends on the parameter cL, 
and thus the closed orbits are possible, but there is no 
potential for which all orbits are closed, since the closing 
depends on cL. 

In the non-relativistic case it is well known that this 
problem has two solutions because the first term in the 
right-hand-side is absent. Then the solutions of eq. (|28p . 
producing closed orbits close to the circular one, are the 
power law potential with V(r)ar@ ~ 2 . For instance the 
natural (3 correspond to the powers -1, 2, 7, 14 ... For 
more general orbits, we can go up to the third order in 
the Fourier series for r/, 

t] = h + h x cos 139 + h 2 cos 2(39 + h 3 cos 3(39 (29) 

and this already produces Bertrand's theorem, stating 
that the closed orbit condition is, 



(3 2 {(3 2 - l){(3 2 -4) = 



(30) 



This includes only the Coulomb and the harmonic oscil- 
lator cases, which indeed have all orbits simply closed for 
a non-relativistic kinetic energy. 

Again, for an ultrarelativistic kinetic energy there is 
no potential with all orbits closed. 



V. CONCLUSION 

Here we study possible quark models for the large de- 
generagy in the meson spectra reported by Bugg [lj. 

We start with a semi-relativistic chiral invariant quark 
model, with relativistic kinetic energy, with negative en- 
ergy components, but with an instantaneous potential. 
For excited states all the different spin-tensor poten- 
tials [13, HH, HH merge and we arrive at a schrodinger- 
like ultrarelativistic linear potential quark model with a 
simple J 2 dependence. Although the chiral degeneracy 
is expected to only dominate the high excitations, this 
schrodinger-like quark model is conveniently simple to 
address the large degeneracy. 

The most ambitious approach to the large degener- 
acy consists in asking for a model with a principal quan- 
tum number. Notice that when the corresponding clas- 
sical orbits are all closed, the Hermann-Jacobi-Laplace- 
Runge-Lenz vector is conserved. In quantum mechan- 
ics this corresponds to a larger symmetry, both angular 
and dynamical, and to a single principal quantum num- 
ber that encompasses both the angular and radial exci- 
tations. Unfortunately we find that this is not possible 
in our ultrarelativistic and instantaneous framework. 

Another possible approach is the one of an approxi- 
mate large degeneracy. The linear and radial excitations 
are so regular in the framework of the ultrarelativistic 
and equal time quark model that approximate patterns 
occur in the spectrum. Notice that in eq. (fTS|) , in Figs. 
[5] and [31 and in Table IIIII an increment of 3 in j is ap- 
proximately equivalent to an increment of 2 in n. For 
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instance the meson with j = 5, n = has a similar mass 
to the one with j = 2, n = 3. Notice that the j = 5 state 
may have / = 4, 5, or 6, while the j ' = 2 state may have 
£ = 1, 1, or 3. we get an approximate but relatively large 
degeneracy of states with different I ranging from 1 to 6. 
To agree with the experiment, we only need to have an 
increment in n much closer to the increment in j. 

Thus, both for the principal quantum number, and 
for an approximate degeneracy, a departure from the ul- 
trarelativistic and equal time quark model is necessary. 
It is plausible that either retardation, string effects, or 
coupled channel effects should improve the theoretical 
model. For instance Morgunov et al. !33[ already showed 
that including the angular momentum of the rotating 
linear confining string changes the slopes of the angular 
Regge trajectories. Moreover, if one could succeed to in- 
clude retardation properly in the confining quark model, 
then one might search again for the possible existence of a 
principal quantum number. And coupled channels might 
also affect differently the angular and radial excitations 



of the spectrum. 

In any case the large degeneracy apparent in the data 
analysis of Bugg, supported by the solid pattern of lin- 
ear angular and radial Regge trajectories, remains quite 
plausible from a theoretical perspective. More experi- 
mental data on the remarkable patterns of excited reso- 
nances is welcome. 
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